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ABSTRACT 

Context. Water is a key tracer of dynamics and chemistry in low-mass star-forming regions, but spectrally resolved observations have so far been 
limited in sensitivity and angular resolution, and only data from the brightest low-mass protostars have been published. 

Aims. The first systematic survey of spectrally resolved water emission in 29 low-mass (L < 40 L©) protostellar objects is presented. The sources 
cover a range of luminosities and evolutionary states. The aim is to characterise the line profiles to distinguish physical components in the beam 
and examine how water emission changes with protostellar evolution. 

Methods. H2O was observed in the ground-state lio-loi transition at 557 GHz (£up/^B~60 K) as single-point observations with the Heterodyne 
Instrument for the Far-Infrared (HIFI) on Herschel in 29 deeply embedded Class and I low-mass protostars. Complementary far-IR and sub-mm 
continuum data (including PACS data from our programme) are used to constrain the spectral energy distribution (SED) of each source. H2O 
intensities are compared to inferred envelope properties, e.g., mass and density, outflow properties and CO 3-2 emission. 

Results. H2O emission is detected in all objects except one (TMCIA). The line profiles are complex and consist of several kinematic components 
tracing diff'erent physical regions in each system. In particular, the profiles are typically dominated by a broad Gaussian emission feature, indicating 
that the bulk of the water emission arises in outflows, not in the quiescent envelope. Several sources show multiple shock components appearing in 
either emission or absorption, thus constraining the internal geometry of the system. Furthermore, the components include inverse P-Cygni profiles 
in seven sources (six Class 0, one Class I) indicative of infalling envelopes, and regular P-Cygni profiles in four sources (three Class I, one Class 0) 
indicative of expanding envelopes. Molecular "bullets" moving at >50 km s~^ with respect to the source are detected in four Class sources; three 
of these sources were not known to harbour bullets previously. In the outflow, the H2O / CO abundance ratio as a function of velocity is nearly 
the same for all line wings, increasing from 10~^ at low velocities (<5 km s~^) to >10~^ at high velocities (> 10 km s~^). The water abundance 
in the outer cold envelope is low, >10~^^. The diff'erent H2O profile components show a clear evolutionary trend: in the younger Class sources 
the emission is dominated by outflow components originating inside an infalling envelope. When large-scale infall diminishes during the Class I 
phase, the outflow weakens and H2O emission all but disappears. 

Key words. Astrochemistry — Stars: formation — ISM: molecules — ISM: jets and outflows 



1. Introduction 

Water is a unique probe of physical and chemical conditions 
in low-mass star-forming regions. In the molecular envelope 
from which the protostar accretes matter, w ater forms almost 
exclusively on the surfa ces of dust grains (iTielens & HagenI 
11982', 'loppolo et al. 2008) and it is the d ominant ice constituent 
dWhittet et al. 1988, Pontop pidan et al.i .2004). As matter flows 
towards the protostar, it heats up a nd water desorbs when the 
dust temperature exceeds 100 K (iFraser et al.l 1200 1) leading 
to a jump in the abundance of several or ders of magnitude 
dCeccarelli et al.|[l996[ iBoonman et al. 2003). At the same time, 
the protostar drives bipolar jets and winds into the envelope. It 
is unknown whether these are molecular or atomic in nature, re- 
gardless of their exact composition, water forms efl&ciently in 
the shocks involved, again enhancing the abundance by several 
orders of magnitude compared to the abundance in the quiescent 



part of the envelope (iKaufman & Neufeldlll996l iFranklin et al.l 



l2008llKristensen et al.ll20TTir 

The water molecule is an asymmetric rotor with many ro- 
tational transitions in the sub-millimetre (sub-mm) and far- 
infrared wavelength ranges. The combination of the abundance 
jumps and the high number of transitions makes water not only 
one of the most important molecular coolants in star-forming 
objects, but also one of the most effective diagnostics for under- 
standing the chemical and dynamical evolution of star forma- 
tion. Altogether, these eflTects make water a much better tracer 
of the various energetic components in a protostar than other 
commonly used molecules such as CO, which do not show large 
abundance variations. 

Water has been observed previously in deeply embedded, 
low-mass protostars, primarily with space -based observatories 
such as ISO, SWAS, Odin and Spitzer (iKessler et al.l I1996L 

' 3- 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



iMelnick et al1l2000LlNordh et al.ll2od3llWemer et al.ll2'QQ4h . ISO 
and Spitzer both observed spectrally unresolved water lines in 
the wavelength range of 5 - 180 pm ( Ceccarelli et al. 19981 
Boosertetal. 2000, Giannini et al. 200 li iLarsson et all I2002L 
Maret et al.. .200Z. .Watson et al.. .2007.) . where many higher- 
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excited water lines (£^up/^B~ few 100 K) were detected. With 
SWAS and Odin it became possible to spectrally resolve 
the ortho-H20 lio-loi ground-state transition at 557 GHz 
(£'up/^B~60 K). The line profiles turned out to consist of a nar- 
row absorption component (Af~ a few km s"^) superposed on 
a sing le outflow co mpon e nt (Al'> 15-2 km s^) (Neufeld et al. 
'2000, Ashbv et al.' '2000", 'Bergin et al." '2003", Hjalmarson^etaL| 
^2003, FrankliiLet al. 2008, B jerkeH et al . 2009). The large beams 
of Odin and SWAS (~2' and 3x4', respectively) were naturally 
more sensitive to large-scale extended emission in the outflow. 

While water emission turned out to be straightforward to de- 
tect in the de eply embedded Class young stellar objects (YSOs; 
iNisini et an i2002), it was not detected in the more evolved Class 
I YSO£]. In particular, neither SWAS nor ISO detected water 
emission from any Class I objects. Two explanations were put 
forward (iNisini et al.ll2002b : either the non-detection is due to 
lower excitation in Class I objects, or it is due to a lower abun- 
dance. If the water is predominantly formed and excited in out- 
flows, the Class I sources with their weaker outflows would go 
undetected, i.e., the water molecules would be less excited. On 
the other hand, the Class I envelopes are usually of lower mass 
and column density compared to their Class siblings, allowing 
for the UV radiation from the accreting star to penetrate deeper 
into the envelope and outflow thereby photodissociating water 
more efficiently. The photodissociation naturally leads to a lower 
water abundance in the Class I objects rendering the emission 
difficult to detect. 

With its much higher sensitivity and spatial resolution 
than any previous instrument capable of observing wa- 
ter, the Heterod yne Instrument for the Far-Infrared (HIFI; 
Ide Graauw et al.l 12010 ) on Herschel makes it possible to de- 
tect water emission from both Class and I sources. Because 
of the improvement in sensitivity and resolution, early HIFI re- 
sults already reveal water line profiles that are much more com- 
plex and show more kinematic components than those observed 
with SWAS and Odin, and than th ose of other species, such as 
CO or SiO (iLefloch et al. l 2010, Kr istensen et al.ll2010bl:[2QTTI) . 
However, this analysis is still limited to only a handful of Class 
sources associated with strong outflow activity, and is therefore 
not representative of how water emission and its profiles evolves 
in young low-mass YSOs, in particular the evolution from Class 
to Class I objects. 

As part of the "Water i n star-forming regions with Herscher 
key programme (WISH; Ivan Dishoeck et al. I l2011h . 29 low- 
mass YSOs were observed with HIFI (see Table [T] for an 
overview). The selected sources span a range of luminosities 
(from -0.8 L© to 40 L©), bolometric temperatures (from ~ 
30 K to ~ 600 K) and envelope masses (-0.01 M© to 10 
Mq), and are all closer th an 450 pc. The Cla ss sources were 
selected from the Hst of lAndre et all (l2000b and the Class I 
sources were se lected from the lists of Tamu ra et al.l ([1991), 
Andre & Mont merl^ (11994 and the S pitzer "Cores to Disks" 
legacy programme (lEvans et al ] l2009l) . Furthermore, these ob- 
jects have been the targets of numerous physical and chemi- 
cal surveys (e.^., [H ogerheiide et al. 1999, Shirlev et al. 2000, 
'J0r2ensenet al. 2002; 2004b|; I2007L lEmprechtinger et al.l 120091 
iKristensen et al. 2010a: b). Some of these objects have also 
been observed with PACS (the Photodetector Array Camera 



^ In the following, Class sources are charactised as YSOs with a 
bolometric temperature, Tboi, < 70 K. The Class sources observed here 
are all Stage sou rces, and the Class I sources are Stage I sources (see 
lEvans etal for a discussion of Tboi as an evolutionary parameter 

and a discussion of the difference between Stage and Class.). 



and Spectrometer) on Herschel, e.g. , HH4 6-IRS and LI 157 
(Ivan Kempen et al.l I201QL INisini et al.l l201Ql) . Eventually HIFI 
and PACS data will exist for all sources in several transitions 
of water, and in other sp ecies such as CO and OH (for a sum- 
mary and overview, see Ivan Dishoeck et al. This paper 
focusses on observations of the ortho-H20 ground-state lio-loi 
transition at 557 GHz in low-mass YSOs, and will thus greatly 
expand on previous SWAS and Odin observations. 

The paper is organised as follows. The observations are de- 
scribed in Sect. [2l Section [3] contains the results, which are all 
discussed in Sect.lH Concluding remarks are to be found in Sect. 

s 

2. Observations and complementary data 

2.1. Herschel-HIFI observations 

A sample of 29 low-mass Class and I YSOs was targeted 
with the Heterodyne Instrument for the Far-Infrared (HIFI; 
"de Graauw et al.' '2010") on the Herschel Space Observatory 
( Pilbratt et al. 2010), see Table [T]for an overview of the sources. 
Observations were performed as single-pointing observations in 
either dual-beam- switch mode with a chop of 3' or position- 
switch mode with a throw large enough to ensure that the ob- 
servations were not contaminated by cloud emission (reference 
positions are provided in Table lA. Il in the Appendix). The lat- 
ter mode was primarily used in crowded regions such as the 
NGC1333 region in Perseus, and the Ophiuchus star-forming re- 
gion. The total telescope time (on -h off + overhead) was 13 min 
for Class sources and 20 min for Class I sources. HIFI was 
tuned to the 557 GHz H2O lio-loi line in the lower sideband of 
Band lb to avoid the strong spur in Band la; this setting does not 
include NH3 Iq-Oq at 572 GHz in the upper sideband. The Wide 
Band Spectrometer and High Resolution Spectrometer backends 
were used for both horizontal and vertical polarizations. 

The beam s ize of HIFI is close to the diffraction limit at 557 
GHz, i.e., 39'' (iRoelfsema et al.ll20T2b . The pointing accuracy is 
better than ~2'', with the caveat that the H- and V-polarization 
receivers are pointing towards slightly different positions in the 
sky (the separation is ~ 6V5). The calibration uncertainty is lim- 
ited by the sideband gain ratio in Band lb, but the overall un- 
certainty is likely < 15%. The velocity calibration is better than 
~100 kHz, or ~0.05 km s"^ at this frequency. 

Data were reduced with HIPE version 6.1 doi feOlOl) and 
exported to CLASS for further analysis. Typical system tem- 
peratures are -85 K. A main-beam efficiency of 0.76 is used 
to convert data from the antenna ten iperature scale (T*) to the 
main-beam temperature scale (T^mb) (IRoelfsema et al .1 120 12). H 
and V polarizations are averaged after individual inspection us- 
ing a simple average. The rms noise in the V-polarization data 
is higher than in the H-polarization data by typically -20%. The 
rms noise is typically 8-13 mK in 0.3 km s"^ bins in the aver- 
aged spectra. 

For the case of dual-beam- switch observations the baseline 
is always very stable and linear baselines are subtracted. The 
position-switched observations show more unstable baselines 
and higher-order polynomials are in some cases necessary to fit 
and subtract the baseline. The spectral ranges are 554.6 - 558.6 
GHz (lower sideband) and 566.6 - 570.6 GHz (upper sideband). 
Detected lines besides H2O include CH3OH (568.566 GHz), CN 
(566.947 GHz) and SO2 (555.666 and 567.593 GHz). 

The reduced spectra are all shown in Figs. [T] (Class 0) and|2] 
(Class I). A number of spectra show absorption against the con- 
tinuum in addition to self- absorption. In the cases where these 
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Fig. 1. Continuum- subtract ed H2O lip-lpi spectra of the observed Class sources. The red dashed line indicates the source velocity 

' et al.l]2QQ2[ 



as derived from C^^O data (|j0rgensen 
unbinned spectra. 



Yildiz et al. in prep.). Some spectra have been rebinned for clarity, insets show 



observations were done as dual-beam- switch (DBS) the data re- 
duction was redone, and the level- 1 reference-position data ex- 
amined for emission. None of the DBS reference spectra shows 
emission, and thus the absorption is believed to be intrinsic to the 
source. Several of the absorption lines have flat bottoms at half 
the measured continuum level, i.e., they are saturated, thus cor- 
roborating that absorption below the continuum level is true. The 
observations of Serpens show a weak absorption feature at ^lsr 
= 1 km s ~\ a feature that also shows up in so me observations o f 
CO 3-2 (Dionat os et"aD l2QlQ) but not afl (Grav es etal]l201Qh . 
This feature probably arises from emission in the reference po- 



sition, and since it is easily identifiable and does not interfere 
with the science goals it is ignored. 



2.2. Complementary data 
2.2.1 . CO emission and Fqo 

In addition to the Herschel data, CO 3-2 data at 345 G Hz have 
been o bt ained at the JCMT and APEX ( van K empen et al. 2006; 
l2009allbl iParise et al.ll2006l Yfldiz et al. in prep.) in the form of 
Nyquist- sampled maps that have subsequently been convolved 
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Fig. 2. Same as Fig.[T]but for Class I sources. 



with a 39'' beam to be directly comparable to the H2O data pre- 
sented here. The CO data are presented in the Appendix (Fig. 
IB. 31 and Table iBJl) and will be further discussed in Yildiz et al. 
(in prep.). 

Carbon monoxide is clearly detected in all sources. The pro- 
files have higher than H2O by more than an order of mag- 
nitude (up to 25 K) and consist typically of a single Gaussian 
component with At; ~ 10 km s"^ with a self-absorption feature 
located at the source velocit}Q. Two sources, LI 44 8 -MM and 
BHR71, also show bullet emission in the CO 3-2 line; this was 



^ Av is used throughout to refer to full width at half maximum and is 
always obtained from Gaussian fitting; Av^^x refers to the full width at 
4-cr rms in 2 km s~^ channels (representative of the full width at zero 



known for L1448-MM (iBachiller et al.lll99Ql) but not reported 
forBHRTl (iParise et al.ll2Q06h . 

Values for the CO outflo w for ce, Fqo, have been assem- 
bled from the literature (Table IBJI in the Appendix). Values of 
Fqo are chosen such that they are all measured using the same 
method and corrections. The corrections consist of taking the 
opacity of the emission in the outflow wings into account and 
correcting for inc lination efl'ects using the correction curves of 
ICabrit & Bertoutl(ll992) . When more than one value exist in the 
literature, the average is used here. Eventually these values will 
be derived in a coherent manner for all sources from our own 



intensity) and is always measured from the data, 
profile function. 



.e., without fitting a 
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Table 1. Source parameters. 



Source 


RA 


Dec 


D 




^bol 


^bol 




References 




(h m s) 


(° ' n 


(pc) 


(km s"^) 




(K) 


(Mo) 




L1448-MM 


03 25 38.9 


+30 44 05.4 


235 


+4.7 


9.0 


46 


3.9 


1,2 


NGC1333-IRAS2A 


03 28 55.6 


+31 14 37.1 


235 


+7.7 


35.7 


50 


5.1 


1,3 


NGC1333-IRAS4A 


03 29 10.5 


+31 13 30.9 


235 


+7.2 


9.1 


33 


5.2 


1,3 


NGC1333-IRAS4B 


03 29 12.0 


+31 13 08.1 


235 


+7.4 


4.4 


28 


3.0 


1,3 


L1527 


04 39 53.9 


+26 03 09.8 


140 


+5.9 


1.9 


44 


0.9 


1,3 


CedllO-IRS4 


11 06 47.0 


-77 22 32.4 


125 


+3.5 


0.8 


56 


0.2 


4,5 


BHR71 


12 01 36.3 


-65 08 53.0 


200 


-4.4 


14.8 


44 


3.1 


6,7 


IRAS 15398' 


15 43 01.3 


-34 09 15.0 


130 


+5.1 


1.6 


52 


0.5 


4,5 


L483 


18 17 29.9 


-04 39 39.5 


200 


+5.2 


10.2 


49 


4.4 


1,8 


Ser SMMl 


18 29 49.8 


+01 15 20.5 


230^ 


+8.5 


30.4 


39 


16.1 


9,10 


Ser SMM4 


18 29 56.6 


+01 13 15.1 


230^ 


+8.0 


1.9 


26 


1.9 


9,10 


Ser SMM3 


18 29 59.2 


+01 14 00.3 


230^ 


+7.6 


5.1 


38 


3.2 


9,10 


L723 


19 17 53.7 


+ 19 12 20.0 


300 


+ 11.2 


3.6 


39 


1.3 


11,12 


B335 


19 37 00.9 


+07 34 09.6 


250 


+8.4 


3.3 


36 


1.2 


1,13 


L1157 


20 39 06.3 


+68 02 15.8 


325 


+2.6 


4.7 


46 


1.5 


1,14 


NGC1333-IRAS3 


03 29 03.8 


+31 16 04.0 


235 


+8.5 


41.8 


149 


8.6 


15,2 


L1489 


04 04 43.0 


+26 18 57.0 


140 


+7.2 


3.8 


200 


0.2 


16,17 


L1551-IRS5 


04 31 34.1 


+ 18 08 05.0 


140 


+7.2 


22.1 


94 


2.3 


11,17 


TMRr 


04 39 13.7 


+25 53 21.0 


140 


+6.3 


3.8 


133 


0.2 


18,17 


TMCIA' 


04 39 34.9 


+25 41 45.0 


140 


+6.6 


2.7 


118 


0.3 


18,17 


TMCl 


04 41 12.4 


+25 46 36.0 


140 


+5.2 


0.9 


101 


0.2 


18,17 


HH46-IRS 


08 25 43.9 


-51 00 36.0 


450 


+5.2 


27.9 


104 


4.4 


19,20 


IRAS 12496 


12 53 17.2 


-77 07 10.6 


178 


+2.3 


35.4 


569 


0.8 


4,21 


GSS30-IRS1 


16 26 21.4 


-24 23 04.0 


125 


+2.8 


13.9 


142 


0.6 


22,23 


Elias 29 


16 27 09.4 


-24 37 19.6 


125 


+5.0 


14.1 


299 


0.3 


24,23 


Oph-IRS63 


16 31 35.6 


-24 01 29.6 


125 


+2.8 


1.0 


327 


0.3 


24,23 


RN091 


16 34 29.3 


-15 47 01.4 


125 


+5.0 


2.6 


340 


0.5 


4,23 


RCrA-IRS5A 


19 01 48.0 


-36 57 21.6 


130 


+5.7 


7.1 


126 


2.0 


4,5 


HHIOO-IRS 


19 01 49.1 


-36 58 16.0 


130 


+5.6 


17.7 


256 


8.1 


4,5 



Notes. Sources above the horizontal line are Class 0, sources below are Class I. Obtained from ground-based C^^O or C^^O observations. 

Measured using Herschel-FACS data from the WISH and DIGIT key programmes (Karska et al. in prep.. Green et al. in prep.). Determined 
from DUSTY modelling of the sources; see Appendix Icl '^^^ Using VLBA observations of a star thought to be associated with the Serpens cluster, 
iDzib et all OplO) derive a distance of 415 pc. The coordinates used in WISH; more accurate SMA coordinates of the sources are 15^43"^02^.2, 
-34°09'06r8 (IRAS15398), 04h39"'13!9, +25°53'2a'6 (TMRl) and 04^39"^35!2, +25Mr44r4 (TMCIA; J0rgensen et al. 2009). 

References. (1) |j0rgensen et al.| ([2007|); (2) iHirota et al.l dMl; (3) iHirota et al.| (|2008|); (4) l Evans et al.l ([20091); (5) iKnude & H0g| fT998|); (6) 
Chen et al. (2008); (7) Seidensticker & Schmidt-Kaler ( 1989); (8) Fuller et al. ( 1995); (9) Hogerheiide et al. ( 1999); (10) Eiroa et al. (2008); (1 1) 
J0rsensen et al. (2002); (12) Goldsmith & Semvak ( 1984); (13) Frerkins et al. ( 1987); (14) Straizvs et al. ( 1992); (15) Bachiller et al. ( 1998) ; (16) 
Brinch et al. (2007); (17) Kenvon et al. (2008); (18) Hoserheiide et al. ( 1998); (19) Velusamv et al. (2007); (20) lNoriega-Crespo et al.l(l2004D : (21) 
iWhittet etal..(J997.) : (22) J0rgensen et al.. (.2009.) : r23) .de Geus et aL (.1989.) : (^24) .Lommen et al.. (.2008.) . 



data (Yildiz et al. in prep.); the literature data are used to obtain 
preliminary results. 

2.2.2. Bolometric luminosity and temperature 

The bolometric luminosity, Lboi, and temperature, Tboi, are mea- 
sured for each source using recent Spitzer-IRAC and MIPS 
fluxes (jEvans et al. 2009) and 2MASS fluxes when available 
(ISkrutskie et al.l l2006). Furthe rmore, Herschel-FAC S fluxes ob- 
tained within the WISH ( van Dishoeck et al. |l201 ih and DIGIT 
(Dust, Ice and Gas In Time; PI: N. Evans) key programmes are 
used, when available, as well as sub-mm and mm fluxes obtained 
from the literature. The data reduction and analysis of the PACS 
data will be presented in two forthcoming papers (Karska et al. 
in prep.; Green et al. in prep.). The newer values of Lboi do not 
change significantly (<20%) for most sources compared to liter- 
ature values. The bolometric temperature, Tboi, is also measured 
for each source using the same data; the availability of the PACS 
data covering the peak of the spectral energy distribution in mos t 
sources has improved the value significantly (lEvans et al 



Using the bolometric temperature to classify sources as either 
Class (rboi<70 K) or Class I (rboi>70 K) implies that Cedl 10- 
IRS4 and IRAS15398 are now Class sources and NGC1333- 
IRAS3 is a Class I source. 



2.2.3. DUSTY modelling 

To quantify the large-scale physical structure of the envelopes 
and to compare the water emission to envelope parameters, each 
source is modelle d using the ID dust radiative transfer code 
DUSTY ( Ivezic & Elitzur 19971) following the prescript ions out- 
lined in lSchoier et al.i (^2002) and'j0rgens en et al.l (|2002|) . Details 
of the procedure are provided in Appendix O The density is as- 
sumed to have a power-law structure, and the temperature is ob- 
tained as a function of position by solving for the dust radiative 
transfer through the (assumed) spherical envelopes, given the lu- 
minosity. The free parameters are fit to the SED and the spatial 
extent of the sub-mm continuum. Furthermore, the mass of each 
envelope is obtained. The mass is calculated at the smallest ra- 
dius where either the dust temperature reaches 10 K or where 
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Table 2. Observed properties of the H2O line profiles. 
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Fig. 3. Gaussian decomposition of the H2O lio-lOl spectra to- 
wards the two sources NGC1333-IRAS2A and LI 157. The indi- 
vidual Gaussian components are shown overlaid on the spectra 
and the residual is plotted below. The saturated absorption and 
inverse P-Cygni profile are not fitted by Gaussians. 

the density drops below 10^ cm"^, in both cases corresponding 
to where the envelope merges with the ambient cloud. 

3. Results 

3.1. H20at557GHz 

The H2O 557 GHz line is detected in all Class sources and all 
but one Class I source, TMCl A, and the full data set is presented 
in Figs. [T] (Class 0) andO (Class I). These observations represent 
the first detection of cold H2O emission from Class I YSOs, i.e., 
emission from the ground-state 557 GHz transition. Integrated 
line intensities are provided in Tabled The dynamical range of 
the integrated intensities is more than two orders of magnitude, 
from -0.1 K km s'^ (Oph-IRS63) to -20 K km s'^ (L1448-MM 
and NGC1333-IRAS4A). On average, the integrated intensity of 
the 1 lo-loi line is almost an order of magnitude weaker in Class 
I sources than in Class sources, explaining why the emission 
was not detected with Odin and SWAS. Peak intensities span 
more than an order of magnitude from -0.03 K to 0.9 K. While 
the Class sources typically have higher peak intensities than 
the Class I sources, there are notable exceptions in both cate- 
gories (Class YSOs with low peak intensity: Cedl 10-IRS4, Ser 
SMM3, L723, L483 and L1527; Class I YSOs with high peak 
intensity: R CrA-IRS5, HHIOO-IRS and GSS30-IRS1). Both R 
CrA-IRS5 and HHIOO-IRS are embedded in the dense photodis- 
sociating region (FDR) created by the Herbig Ae star, R CrA 



Source 


rms^ 




C /-T-f total 1 h 






(mK) 


(K) 


(Kkms-^) 


(kms-^) 


L1448-MM 


12 


0.32 


18.7 


138 


NGC1333-IRAS2A 


12 


0.31 


5.2 


58 


NGC1333-IRAS4A 


12 


0.90 


17.4 


84 


NGC1333-IRAS4B 


12 


0.69 


10.6 


54 


L1527 


10 


0.14 


1.2 


22 


CedllO-IRS4 


8 


0.03 


0.3 


11 


BHR71 


12 


0.44 


10.2 


131 


IRAS15398 


8 


0.21 


2.0 


28 


L483 


12 


0.14 


3.2 


22 


Ser SMMl 


12 


0.58 


10.8 


69 


Ser SMM4 


12 


0.47 


8.8 


80 


Ser SMM3 


12 


0.13 


3.6 


80 


L723 


11 


0.13 


0.7 


15 


B335 


13 


0.09 


1.1 


26 


LI 157 


13 


0.17 


3.0 


67 


NGC1333-IRAS3 


11 


0.20 


5.2 


54 


L1489 


11 


0.08 


0.8 


41 


L1551-IRS5 


9 


0.12 


0.7 


19 


TMRl 


8 


0.04 


0.3 


9 


TMCl A 


8 


0.04 


<0.1 




TMCl 


11 


0.05 


0.5 


19 


HH46-IRS 


8 


0.15 


1.2 


24 


IRAS12496 


8 


0.08 


0.6 


24 


GSS30-IRS1^ 


10 


0.72 


4.4 


35 


Elias29 


10 


0.14 


2.0 


28 


Oph-IRS63 


10 


0.02 


0.1 


5 


RN091 


10 


0.16 


0.5 


13 


RCrA-IRS5A 


17 


0.50 


3.8 


34 


HHIOO-IRS 


16 


0.63 


5.3 


28 



Notes. Measured in 0.3 km bins. Integrated over the entire line, 
including bullet emission. Spectra smoothed to 2 km s~^ resolution 
and line width measured at the 4-cr level. Possible contamination 
from the VLA1623 outflow. 




Broad Medium Bullet 



Fig. 4. Histogram over detected features in Class (red) and 
Class I (blue) sources; the features include the broad and 
medium components, bullets and inverse and regular P-Cygni 
profiles (IPC and RPC, respectively). 



(Neuhauser & Forbrichll2008 V which may be causing the higher 
H2O intensity in these two sources. Indeed, both profiles show 
a strong narrow component emission component (Av ~ 3-5 km 



). 



The main diff'erence between the Class and Class I YSOs 
leading to the dichotomy in integrated intensity is the linewidth. 
The linewidth at zero intensity, here defined as the 4-cr level, is 
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Fig. 5. Distribution of Av of the different Gaussian components. 
Class sources are marked in red, and Class I sources in blue. 
The horizontal lines show the range of velocities adopted for the 
narrow, medium and broad components. 



also listed in Table [2l and the range extends over more than one 
order of magnitude. The maximum linewidth, Af;niax, includes 
contributions from bullets. The line pr ofiles are very different 
from envelope-only mod el predictions (IPoelman & van der TakI 
l2QQ7l Ivan Kempen et al.| [2QQ8). indicating that very little emis- 
sion is coming directly from the envelope, but that the bulk arises 
from outflows. In general, very little water emission originates 
in col d, quiescent regions (e.g ., lBergin et"aDl2QlQilCaselli et~an 
l201QL lHogerheii de et al. 11201 ih . 

The most prominent features of the line profiles are their 
complexity and width, with no two profiles being identical. 
Nevertheless, it is possible to identify a set of common traits 
based on Gaussian decomposition of the line profiles (see Fig. [3] 
for examples and Tableland Fig. [4] for an overview): 

- a broad component with Av > 20 km s"^; this component is 
present in most Class sources and several Class I sources 
and when it is detected it is always in emission. The center 
of this component is sometimes shifted by up to 10 km s"^ 
with respect to the source velocity, as determined from C^^O 
lines (e.g., NGC1333-IRAS2A); 

- a medium component with a Av of 5-20 km s"^; in 
NGC1333-IRAS3 and Ser SMM3 this component is seen 
in absorption, otherwise it is in emission. The component 
is sometimes blue-shifted by up to 10 km s"^ with respect to 
the source velocity (e.g., NGC1333-IRAS4A); 

- a narrow component with a Av of ~2-3 km s"^ ; this compo- 
nent is mostly seen in absorption, but not always. When the 
component is in absorption, it is often saturated, i.e., all of 
the continuum photons are absorbed as well. A few sources 
show multiple absorption components at velocities shifted by 
up to a few km s"^ from that of the source (e.g., the Serpens 
sources); 

- extremely-high- velocity (EHV) components or "bullets" ap- 
pearing as near-Gaussian profiles separated from the broad 
or medium components as well as from the source velocity 
by several tens of km s"^; only three Class sources show 
these components in the 557 GHz line (LI 44 8 -MM, BHR71 
and LI 157); 

- regular or inverse P-Cygni-type profiles where narrow emis- 
sion and absorption are seen next to each other and always 
centred at the source velocity. 



The results of the Gaussian decomposition are provided in Table 
IA.2I in the Appendix, and all features are discussed in Sect. 
14.11 Many of these features have never been identified in these 
sources in other species before, such as the bullet emission from 
BHR71 and LI 157. The broad component is observed in sev- 
eral low-mass YSOs, but it never appears so prominently in, 
e.g., CO rotational lines. For the sources showing both medium 
and broad components, the decomposition was checked against 
data for higher- excited lines where the a bsorption feature is not 
interfering (see iKristensen et ani2010aL for examples). Inverse 
P-Cygni profiles have so far only been detected in one low- 
mass YSO, N GC1333-IRAS4A, in H2CQ, ^^CQ, HCN and 
H2Q emissio n Francesco et al ]|2001L togensen et al.ll2007l 
lAttard et al.1 12009. Kristen sen et al.ll2010bh . Our H9O data in- 
crease the sample by almost an order of magnitude. 

Figure [5] shows a histogram of Av of each component. The 
broad H2O component is present in almost all Class sources 
(except CedllO-IRS4, IRAS15398, L483 and L723) but less 
than half of the Class I sources (NGC1333-IRAS3, L1489, 
L1551-IRS5, TMCl, HH46-IRS and IRAS12496). The medium 
component is found in ~2/3 of all the sources, independent of 
evolutionary stage. The narrow component, however, is present 
in nearly all sources. With the exception of two sources (B335 
and HHIOO-IRS), narrow emission always appears together with 
absorption either in a regular or inverse P-Cygni profile. Two 
sources show absorption broader than 5 km s"^ (NGC1333- 
IRAS3 and Ser SMM3). NGC1333-IRAS3 shares a reference 
position with the three other NGC1333 sources, and since the 
absorption is not seen in any of them, it likely originates in the 
source. Similarly, Ser SMM3 shares a reference position with 
SMMl and SMM4, where the absorption is not seen. These ref- 
erence positions were carefully chosen to be free of ^^CO emis- 
sion and thus should be free of H2O emission given the diffi culty 
in detecting any H2O emission from cold clouds (.Caselli et al.l 
I2OIOI) . 



3.2. Correlations between H2O and outflow/envelope 
parameters 

To quantify the origin of H2O emission further and to use the 
full sample of 29 sources, the total H2O integrated line intensity 
(including bullet contributions) is compared to different enve- 
lope and outflow parameters. The parameters are the envelope 
mass (Menv), the bolometric luminosity and temperature (Lboi 
and Tboi), the envelope density at a distance of 1000 AU from the 
source (^(1000 AU)), integrated CO 3-2 emission, linewidth of 
the 557 GHz line (Af;niax), CO outflow force ( Fco) and the alte r- 
native evolutionary parameter L^-^^/M^^y ( Sara ceno et al.lll99"6l) . 

For these comparisons, the observed line intensities are 
scaled by the square of the source distance to a common dis- 
tance of 200 pc, the average source distance. For the compari- 
son with CO 3-2 emission, neither H2O nor CO 3-2 emission 
is scaled. The comparisons are quantitatively examined using 
Pearson's correlation coefficient, r. For values of r close to 1, the 
correlation is good, for values of r close to it is non-existing, 
and for values close to -1 there is a strong anti-correlation. In 
the following, a strong correlation is defined as |r|>0.7 and weak 
correlations are defined to have 0.6<|r| <0.7; these values cor - 
respond to >3.5(T and ~3(T, respectively (iMarseille et al.ll20TQl) . 
The results of all the correlations are shown in Table |4]and sam- 
ple correlations are shown in Fig. [S] where the sources with spe- 
cial features, such as inverse P-Cygni profiles, are highlighted 
separately. 
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Fig. 6. Correlation plots showing the integrated H2O lio-loi 557 GHz line intensity scaled to a distance of 200 pc as a function of 
various envelope parameters; the H2O intensity is not scaled when comparing to CO 3-2 emission. Class sources are marked in 
red, Class I sources in blue; sources with bullet emission are marked with crosses, inverse P-Cygni profiles with squares and regular 
P-Cygni profiles with triangles; the upper limit for TMCIA is marked with an arrow. Strong correlations (|r| > 0.7) are marked with 
full lines, tentative correlations (0.6 < |r| < 0.7) with dashed lines. The parameters are (from left to right, top to bottom): CO 3-2 
emission, Lboi, T'boi, ^env, the evolutionary parameter L^^^^/Menv, ^co, ^(1000 AU) and At;niax- 
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Table 3. Observed features in the H2O lio-loi 557 GHz line profiles^ . 



Source 


Broad^ 


Medium^ 


Narrow^ 


Bullets Saturated abs. 


Multiple abs. 


Inverse P-Cygni 


P-Cygni 


L1448-MM 


E 




A 


X 








NGC1333-IRAS2A 


E 


E 


A 


X 








NGC1333-IRAS4A 


E 


E 


A+E 


X 




X 




NGC1333-IRAS4B 


E 


E 


A 


X^ X 








L1527 


E 




A+E 






X 




CedllO-IRS4 




E 


A 


X 








BHR71 


E 




A+E 


X 




X 




IRAS 15398 




E 


A+E 


X 




X 




L483 




E 


A 


X 








Ser SMMl 


E 


E 


A 


X 


X 






Ser SMM4 


E 




A+E 


X 


X 


X 




Ser SMM3 


E 


A 


A 


X 


X 






L723 




E 


A+E 


X 






X 


B335 


E 




E 










L1157 


E 




A+E 


X 




X 




NGC1333-IRAS3 


E 


A 


A 




X 






L1489 


E 




A 










L1551-IRS5 


E 




A+E 


X 






X 


TMRl 




E 


A 










TMCIA 
















TMCl 


E 














HH46-IRS 


E 




A+E 








X 


IRAS 12496 


E 




A+E 


X 






X 


GSS30-IRS1 




E 


A+E 


X 




X 




Elias 29 




E 


A+E 


X 








Oph-IRS63 




E 












RN091 




E 












RCrA-IRS5A 




E 


A+E 










HHIOO-IRS 




E 


E 











Notes. Detected components marked with E (emission) or A (absorption). Detected features are marked with X. At; > 20 km s ^ 5 < Au 
< 20 km s~^ Af < 5 km s~^ The bullet emission is only detected in the excited 2o2-lii line at 988 GHz. 



Table 4. Pearson's correlation coefficient for H2O and CO emission, envelope and outflow parameters. 



H^O CO TZ\ 7W M Z?^7m F^o KIOOO AU) M 



H20 


0.64 0.56 


-0.60 


0.75 


-0.61 


0.50 


0.71 


0.90 


CO 


0.56 


-0.44 


0.71 


-0.56 


0.46 


0.63 


0.66 


^bol 




0.03 


0.59 


-0.18 


0.28 


0.19 


0.46 


T'bol 






-0.45 


0.56 


-0.70 


-0.66 


-0.49 


M 








-0.90 


0.43 


0.70 


0.58 












-0.38 


-0.75 


-0.46 


Fco 












0.65 


0.51 


^(1000 AU) 














0.61 


A^max 

















H2O emission at 557 GHz correlates the strongest with Af;niax 
(r = 0.90) followed by envelope mass (r = 0.75), and ^(1000 AU) 
(r = 0.71). Emission correlates weakly with CO 3-2 emission 
(r = 0.64), L^;fj/Menv (r = -0.61) and T^oi (r = -0.60), but not 
with Lboi and Fqo- The strong correlation between the H2O inte- 
grated intensity and Af;niax is not a surprise, considering that most 
of the emission of each profile is located in a single medium or 
broad Gaussian component, i.e., the integrated area of the pro- 
file is directly proportional to the peak intensity (Tpeak) times 
the width. To demonstrate that the profiles are dominated by a 
single Gaussian component, the width is measured at 25% of 
^peak; sources with rpeak/4 < 4cr are not included. For a Gaussian 
function, the area is equal to 0.53 x T^peak ^ Al'(25%). Figure 
[7] shows the correlation between J Tub / (rpeak x Af;(25%)) 
and Ai;(25%). The correlation itself is not strong, the scatter be- 



ing dominated by profiles with a width of < 20 km s"^ (Fig. [7]). 
However, with the exception of two sources, all points have val- 
ues close to the expected, 0.53, and we conclude that the profiles 
are indeed dominated by a single Gaussian component. 



The correlations are displayed in Fig. [8] and discussed in 
Sect, m To further confirm that the correlations are real, the 
partial correlation facto rs are calculated for a ll combinations of 
three parameters (e.g., iMarseille et al.ll2QlQl) . At the 2cr level, 
all strong correlations are confirmed except that between den- 
sity and envelope mass, i.e., there is a 95% likelihood that the 
correlations are real. 
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Fig. 7. Correlation between the width measured at the 25% level 
of the peak intensity, Tpeak and the integrated intensity ( J Tmb 
df;) divided by Tpeak- Symbols are as in Fig. (6] The thick dashed 
line is for a proportionality constant of 0.53, corresponding to a 
Gaussian curve. 



Fig. 9. Gaussian decomposition of the source Ser-SMM3. This 
source shows a broad emission feature (At; ~50 km s"^). A 
medium absorption feature (Ai; ~15 km s"^) is superposed on the 
emission along with two narrow absorption features at i^lsr~1 
km s"^ and ~8 km s"^ The former of the narrow absorption fea- 
tures is likely due to emission in the reference position (see text) 
and the latter is intrinsic to the source and is saturated. 



H2O 




Fig. 8. Schematic diagram showing which parameters correlate 
with each other. Only strong correlations with |r| > 0.70, for 
which the correlation is confirmed, are shown. 



4. Discussion 

4. t. H2O and dynamics: profile components 
4.1 .1 . Outflow components 

The broad and medium line components have widths exceed- 
ing the line width typically seen originating from the quiescent 
envelope (~ 1 km s"^; e.g.. *j0rge nsen et aP 1200 2) and are thus 
likely related to the protostellar outflow. The justification for dis- 
tinguishing between two components comes from NGC1333- 
IRAS2A, NGC1333-IRAS3, NGC1333-IRAS4A, Ser SMMl 
and Ser SMM3, where both components are seen, and they 
are off'set from each other by up to 10 km s"^. The ofl'set 
results in ver y asymmetric line profiles, as also illustrated in 
iKristensen et al.^ (2010b) (see Fig. [3] and O. The statistics pre- 
sented in Fig. [5] show that the broad component is more closely 
associated to the earlier evolutionary stage, whereas the medium 
component is equally present at both stages. 



Viewing angle. The viewing angle may play a role in whether 
a source has a broad component associated with it, or both a 
broad and medium component, but because of the relatively 
large sample size such an eff'ect should be statistically small. 
For example, the two Class sources L483 and LI 527 both 
have outflows moving close to the plane of the sky and none 
of them show broad emission components. Yet, the outflow in 
NGC1333-IRAS2A is likely moving close to the plane of the 
sky, but the profile shows a broad component. Because the Class 
I sources show weaker emission it is also conceivable that broad 
emission components are hiding in the noise in these more 
evolved sources. For the weak Class I sources (TMRl, Oph- 
IRS63, RN091) this may be the case, but the stronger Class 
I sources (GSS30-IRS1, EHas29, RCrA-IRS5 and HHIOO-IRS) 
should have detectable broad components associated with them 
if they exist. Therefore, the disappearance of the broad com- 
ponent at later evolutionary stages is likely real, whereas the 
medium component remains through the Class I stage, but de- 
creases in intensity. The disappearance implies an evolutionary 
trend in the broad component as it weakens in both width and 
intensity from the Class phase to the Class I phase. 

Broad or medium. In sources where only the broad or the 
medium component is observed, the profile typically corre- 
sponds to what has previously been observed in low-/ CO 
emission from ground-based telescopes. This component is spa- 
tially extended over > 1000 AU and associated with the large- 
scal e outflow. A m ap of H2O emission in LI 157 from RACS 
(Nisi ni et aP 12010 ) shows that this component is physically 
closer to the region which is currently being shocked (as traced 
by H2 emission) than to the swept-up gas in the molecular out- 
flow (as traced by low-/ CO emission). With the data presented 
here, it is not possible to infer whether H2O emission is always 
associated with shocked gas rather than with swept-up gas. No 
matter the exact excitation conditions, when only the broad or 
medium component is present it traces the outflow. 
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Broad and medium. In the sources where both a broad and 
medium component is present (NGC1333-IRAS2A, NGC1333- 
IRAS3, NGC1333-IRAS4A, NGC1333-IRAS4B, Ser SMMl 
and Ser SMM3), the medium component is either centred at the 
source velocity (offset < 1 km s-^NGC1333-IRAS2A) or blue- 
shifted by 2-10 km s"^. The medium component is not seen in 
CO 3-2 or any other species, except for the case of NGC1333- 
IRAS2A. Furthermore, the medium component is seen in ab- 
sorption against the broad component in NGC1333-IRAS3 and 
Ser SMM3 (see Fig. [9] for the example of Ser SMM3). The ab- 
sorption implies that the medium component is located in front 
of the blue- shifted part of the broad component, and that the 
medium component is less extended than the broad (if it were 
more extended, it would also show emission; see Fig. [TOlfor a 
cartoon illustration). The broad component is not observed in 
absorption, neither against the continuum nor the medium com- 
ponent. 

If there is a corresponding red- shifted medium component 
(not seen in any of the objects) it must be located behind the 
red-shifted part of the broad component, and be less extended 
than the broad component. Furthermore, the broad component 
must be optically thick over the extent of the red- shifted medium 
component, if it exists, or the corresponding red-shifted medium 
component would have been detected in emission. The medium 
component must be located inside the envelope; the narrow ab- 
sorption components are saturated in Ser SMM3 and NGC1333- 
IRAS3 showing that the absorbing envelope is between us and 
the medium absorption component. These geometrical consid- 
erations along with the non-detection of the component in CO 
point to an origin in the inner part of the envelope on scales < 
1000 AU. 

Viewing angle appears not to play a role in determining 
whether a source shows both a broad and medium component. 
The three NGC1333 sources IRAS2A, IRAS4A and IRAS4B 
have outflows moving close to the plane of the sky, close to 
an angle of 45° and close to the line-of-sight, respectively 
(|j0rgensen et al.ll2Q04aL Yildiz et al. in press). 



Association with masers. A Gaussian peak off'set from the 
source velocity is a signature of a Jump-type shock, as opposed 
to a Continuous-type shock which s hows up as a trian gular pro- 
file peaking at the source velocity (Hollenbach 1997). This is 
not the first time H2O emission is linked with J-type shocks 
on small spatial scales in protostellar envelopes; H2O masers at 
22 GHz are also thought t o be ass ociated with J-type shocks in 
dense gas (see review by lElitzu3 il992). Many low- mass pro- 
tosta rs have H2O masers associated with them (Furuy a et al.l 
120031 and VLBI observations show them to originate typically 
within the inner few 100 AU of the source. The sources show- 
ing both broad and mediu m componen ts all have H2O masers 
associated with them (e.g. jFuruya et al . 2003, Moscadelli et al.l 
^06). Most of the masers are variable in terms of velocity over 
time-scales of years or less. It is therefore no surprise that maser 
spots are found matching the observed velocity of the medium 
components seen here. It is interesting to note that in some of 
t he sources w ith little maser variability (e.g., NGC1333-IRAS3; 
[Furuya et al.ll2003l) . the velocity of the maser spot coincides with 
the velocity of the medium component (4.6 km s"^ versus 4.9 km 
s~^). The same is seen in Ser S MMl in the 183 GHz H2O maser 
line dvan Kem pen et al."2009cV where the velocity of their maser 
spot 2 is 7.4 km s"\ identical to the velocity of the medium com- 
ponent in Ser SMMl. It is important to note here that we do not 
propose that the maser causes the medium component. However, 



the phenomenon producing the masers likely has a thermal com- 
ponent associated with it, which produces the medium compo- 
nent. This is similar to what ha s previously be en ob served in 
high-mass star-forming regions (iMarseille et al.l 120101) . Indeed, 
the physi cal cond itions for maser emission, ^ > 10^ cm"^ and T 
~ 400 K (I Elitzuii [T992). are such that a thermal component must 
be present. The only limiting factor is the column density of the 
thermal component: if the beam-averaged column density is low, 
the component is not detected. 

Five sources have H2O maser activity associated with them, 
but do not show both broad and medium components (Ser 
SMM4, L483, L723, LI 157 and GSS30-IRS1). The water 
masers a re not constant ove r the observed epochs in these five 
sources (iFuruya et al.l l2003l) . but seem to switch on and off*. 
When the masers are on in these five sources, the intensity is 
weaker than in the sources showing both broad and medium 
components. Furthermore, L483, L723 and LI 157 are weak in 
557 GHz emission, and it cannot be ruled out that the non- 
detection of the medium component is a question of S IN. 

While the exact origin of the masers and the J-type shock 
causing the medium component is, is unknown, however there 
are several candidate mechanisms. If the protostellar wind di- 
rectly impacts on the cavity wall, this could drive a J-type shock 
into the cavity wall, a mechanis m which has bee n proposed for 
masers in high-mass protostars (God di et al.l2006h . Another pos- 
sibility is that the masers origi nate in the jet itself, or at least that 
they are connected to the jet dFuruya et al.ll2003l) . which is al- 
most certainly true for the high- velocity masers (f;>20 km s"^). 
Because the medium component is located physically between 
the inner outflow close to the base and the outer absorbing enve- 
lope, these observations suggest that the lower- velocity masers 
are more directly linked to the outflow cavity walls rather than 
the jet itself, if the physical mechanism behind the masers and 
the medium component is the same. 

4.1 .2. Molecular bullets 

Molecular "bullets" or extremely high- velocity (EHV) gas 
has been detected in four Class sources, LI 448, BHR71, 
NGC1333-IRAS4B and LI 157, with the bullet emission in 
NGC1333-IRAS4B only seen in the higher- excited H2O 2o2- 
lii line at 988 GHz (iKristensen et al.ll2010bl) . The three previ- 
ously unknown bullet features highlight H2O as a good dynam- 
ical probe, ther eby co mplementing p revious SiO observations 
of bullets (e g..lGirart& Acord 200 11 iNisini etal.ll2002l I2007L 
ICabrit et al.1 l2007l) . IKristensen et al.l (l201ll) analysed the H2O 
bullet emission in LI 448, and came to the conclusion that the 
bullets are rich in H2O and H2, the former requiring the pres- 
ence of the latter. In particular, H2O abundances of lO'^-lO""^ 
are derived with respect to H2. 

High spatial resolution images of, e.g., CO and SiO, show 
that the bullets consist of sub-arcsecond-sized clumps ty pjca lly 
separated by 1-2'' in the plane of the sky (e.g., iMaury et aP 
2010, for the case of LI 448). The bullets are located along the 
jet axis, wher e they are likely internal working surfaces along 
the jet (e.g., Santiago-Garci a et al.ll2009l) . 

The bullets appear in more sources than were previously 
known to harbour them, indicating that the phenomenon is more 
common than previously known. Whether all protostars show 
bullet emission, whether it is a specific phase that protostars go 
through or whether it is limited to specific sources will need 
to be tested by observing more sources at higher S IN and in 
other species. To date, bullets have not been detected in Class 
I sources, and it is not known whether this is due to an ab- 
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Fig. 10. Cartoon showing the different dynam- 
ical and energetic components in protostars di- 
rectly traced by water emission. Each panel 
shows a cartoon of a protostellar system with 
a sample spectrum and a number of specific ex- 
amples of objects in this category. Note, the ve- 
locity scale of the inserted spectra change from 
the top left panel to bottom right. Inverse P- 
Cygni and regular P-Cygni profiles have been 
abbreviated IPC and RPC, respectively. The la- 
bels Class and I refer to the Class where this 
type of profile is typically found. 



sence of the mechanism creating the bullets, or if it is a lack of 
S IN. For the two sources showing the brightest bullet emission, 
L1448-MM and BHR71, the bullet emission contributes -25- 
30% of the total integrated intensity of the 557 GHz emission 
(for LI 157 the bullet contribution is ~15%). Furthermore, the 
bullets are only seen in Class sources with a broad component. 
Of the Class I sources with a broad component, only NGC1333- 
IRAS3 is bright enough that bullets would be detected, were they 
present. Therefore it is not possible to conclude whether the phe- 
nomenon is limited to Class sources; deeper integrations on 
Class I sources are required to determine whether it is particular 
to Class sources. 



4.1 .3. Absorption features 

The narrow component (At; < 5 km s"^) is seen in absorption 
in all but one Class source and a few Class I sources (B335, 
TMCl, Oph-IRS63, RN091 and HHIOO-IRS). This component 
is interpreted as being caused by absorption by the outer enve- 
lope and ambient cloud, consistent with previous interpretations 
dKristensen et al.|[2010bh and originates in the same region a s 
low-J C^^O and HCO+ emission (e.g., |j0rgensen et al.(l2002h . 
Because the absorption is seen against both the outflow and con- 
tinuum emission, the absorbing layer must be located in front of 
both the emitting layers, i.e., the outflows are also embedded. 
If the absorption layer were located behind the outflow, flat- 



bottomed absorption features not extending down to the zero- 
level should have been observed; this is not the case. 

The non-LT E radiative transfer code RADEX 

(Ivan der Tak et al.l l2Q07l) is used to estimate the H2O col- 
umn density required to cause the saturated absorption. For 
physical conditions such as those found in the outer envelope 
(^H=10^ cm"^, r=15 K, Af;=l km s"^) an 0-H2O column 
density of > 10^^ cm"^ is required to reproduce the saturated 
absorption, i.e., for those conditions, the line optical depth, r, is 
more than 3. The result is not sensitive to the precise physical 
conditions since the bulk of the population is in the loi ground 
state. Flat-bottomed profiles typically require r > 7, which 
is obtained for a column density of -3-4x10^^ cm"^. Thus, 
deep non-saturated profiles are found in a very narrow range of 
column densities. Limits obtained from deep H2^0 observations 
will be used to further constrain the column density of H2O in 
the outer envelope. 

The H2 column density calculated from spherical en- 
velope models ranges from a few times 10^^ cm"^ (e.g., 
L483, L1489, TMCl, IRAS12496) to -10^4 cm'^ (NGC1333- 
IRAS4A, NGC1333-IRAS4B) in a pencil beam towards these 
sources (see Appendix O). Therefore, the 0-H2O abundances 
range from > 10"^^ for the highly embedded sources, to >10"^ 
for the less embedded sources. To determine the abundance of 
the absorbing layer more accurately, envelope models are re- 
quired taking the range of densities and temperatures found in 
the envelope into account. NGC1333-IRAS2A was modelled in 
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some detail (iKristensen et al.|[2QlQah . where the abundance is of 
the order of 10"^. The abundance is higher than what is found for 
other Quiescent regi ons such as prest ell ar core s or protoplanetary 
disks (iCaseUi et al.ll201Ql TBergin et alJl201Ql iHogerheiide et alJ 
l2QTTh . Because these values are all lower limits, there is no rea- 
son to believe that there is a real trend with evolution or col- 
umn density. There is, however, a trend for more Class sources 
showing saturated absorption (10/15) than Class I sources (4/14). 
If the abundance in the outer envelope is the same in Class 
and I sources, then this difference can arise from the Class 
envelopes having a higher H2 column density than the Class I 
sources. In this scenario, the absorption is less likely to be satu- 
rated in the more dilute Class I envelopes. 

4.1.4. Inverse P-Cygni profiles 

Inverse P-Cygni profiles are detected towards seven sources: 
NGC1333-IRAS4A, Ser SMM4, LI 157, IRAS15398, L1527, 
BHR71 (all Class 0) and GSS30-IRS1 (Class I). An inverse 
P-Cygni profile is characterised by a narrow absorption fea- 
ture red-shifted with respect to a narrow emission feature and 
is typically caused by a velocity and excitation gradient in gas 
falling in towards a central object, in this case the protostar. 
The profiles are located on top of the broader outflow compo- 
nent, and the absorption takes place against both the outflow 
component and the continuum. This implies that the absorb- 
ing layer is located between the outflow component and the 
edge of the envelope. Inverse P-Cygni profiles have previously 
been detected for NGC1333 -IRAS4A, in both interferometric 
and single-dish da ta ( Di Fran cesco et a l.' '2001', 'J0rg ensen et al.l 
I2OO7L lAttard et al.ll2Q09. Kristensen et al. 2010b) . These obser- 
vations constitute the first detections of seven inverse P-Cygni 
profiles in a 40'' beam of the envelopes surrounding low-mass 
protostars. As noted above (Sect. O, none of the observations 
show signs of contamination from the reference position or fore- 
ground material. 

To quantify the infall in a simplistic manner, 
iDi Francesco et al.l (l200ll) used a toy model consisting of 
two slabs moving towards each other with a continuum source 
between them by modifying the model of Myers et al.l 11996). 
The two slabs are considered to be isothermal and of uniform 
density. There are two limitations to the model: the simplicity 
and the large number of free parameters. First, an inf ailing 
envelope has temperature, velocity and density gradients along 
the line of sight which are not taken into account. Second, the 
large number of free parameters for such a toy model includes 
the excitation temperature of the front layer, the rear layer 
and the continuum source (Tf, T^, T^), the filling factor of the 
absorbing source with respect to the telescope beam (O), the 
peak optical depth of the front and rear layers (tq), the turbulent 
velocity dispersion in each layer (Af;tur) and the infall velocity 
(f;in). Furthermore, the absorption also takes place against the 
outflow component in Ser SMM4, NGC1333-IRAS4A and 
GSS30-IRS1, where the absorption is saturated. Absorption 
against multipl e layers is taken in to account in the modified 
4-slab model of lAttard et al.l (l2009l) . Introducing four slabs with 
different velocities and excitation temperatures increases the 
number of free parameters significantly, making any solution 
even more degenerate. Instead of adopting more parameters to 
reproduce the entire profile, only the outflow- subtracted inverse 
P-Cygni profiles are analysed in the following. The outflow 
contribution only plays a part in the value of O, which is larger 
than if just the continuum emission is considered, and the value 
of Tc, which includes both a contribution from the source and 
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Fig. 11. Inverse (top) and regular (bottom) P-Cygni profiles de- 
tected in the eleven sources (red) with the infall or expansion 
model shown on top (black). The broad outflow component has 
been subtracted from each spectrum. Spectra have been centred 
at km s"^ 



outflow. The main focus will be on the dynamics rather than the 
excitation. 

Any infall velocity found here should be considere d as 
a "characteristic" infall velocity as also noted in Myers e t al.l 
(1996). To estimate the excitation conditions, H2O is assumed to 
be sub-thermally excited because the density is likely lower than 
the critical density for this transit ion, -10^ cm"^ in the optically 
thin limit (iDubemet et al.ll2QQ6l) . For densities in the range of 
10^ - 10^ cm~^, the excitation temperature in the front layer is 
3-4 K for a column density of 10^^ cm"^ (Sect. 14. 13]) and a ki- 
netic temperature of 10-20 K. Under the same assumptions, the 
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excitation temperature of the rear layer is c loser to the kinetic 
temperature (see Fig. 1 of iMyers et al.lll996[ for details), which 
in this case is expected to be 10-20 K for the outer parts of the 
envelope. The excitation temperature of the continuum layer is 
assumed fixed at 15 K; this parameter is degenerate with and 
To, and it is possible to find other solutions if the continuum tem- 
perature is changed to, say, 10 or 20 K. However, since the goal 
is not to find the best-fit solution for the excitation conditions, 
this is ignored at present. 

The shape of the profile is dominated by f^inf/Af;turb, and it 
is therefore possible to estimate the infall velocity of the two 
slabs along with the turbulent width, Af;turb • The fitting is done 
by eye, as the goal is to get a qualitative estimate of the infall 
velocity. The numbers are given in Table |5] along with the other 
parameters used to describe the profiles, and the resulting pro- 
files are compared to observations in Fig. [TT] Characteristic in- 
fall velocities range from 0.1 km s"^ (LI 157) to 1.2 km s"^ (Ser 
SMM4) for Af;turb ranging from 0.2 km s"^ (BHR71) to 1.0 km 
s"^ (Ser SMM4, L1527). The turbulent velocities of Af;(C^^O 
3-2) agree to within 0.1-0.3 km s"^ with the values found here 
dlirgensen et al. 2002, ; Yildiz et al. in prep.). The uncertainty on 
both the infall velocity and the velocity dispersion is estimated 
to be 0.2 km s~^ by comparing difi'ere nt m odel profiles to the 
observations. iDi Francesco et al.l (l200lh and lAttard et al.l (l2009l) 
both modelled the infall in NGC1333-IRAS4A using the same 
method but for diff'erent molecules (H2CO and HCN), and both 
found a typical infall velocity of 0.6 km s"^ in good agreement 
with 0.5 km s"^ inferred here. 

If the envelopes are in free fall thro ughout, th e velocity at 
any point through the envelope is t; = ^f2GM^^Jr, where Menc 
is the ma ss enclosed within radius r. In the classic inside-out Shu 
collapse (IShull977h . the infall wave travels through the envelope 
at the sound speed, and only the inner parts of the envelope are 
free-falling towards the centre at early times. If the inferred en- 
closed mass is less than the total envelope mass obtained from 
DUSTY modelling, the infall wave has likely not reached the 
outer edges of the envelope, and the assumption that the entire 
envelope is in free fall is no longer valid. The infall velocity at 
the outer edge of the envelope will therefore range between 
km s"^ and the free-fall velocity at the outer radius. Moreover, if 
the envelope is rotating, the infall velocity may be systematically 
lower than the free-f all velocity, depending on the viewing angle 
(iTerebev et al .""1984'). 

The characteristic distance of the infall zone for the en- 
velopes that may be undergoing free-fall collapse can be esti- 
mated. For infall onto a 0.5 M© central object the velocity is 
6 km s"^ at 25 AU. This stellar mass is a typical upper limit 
given that Class sources are thought to have accreted less than 
half of their final mass ( Andre et al. 2000), and the peak of the 
initial mass function is at 0.3-0.5 re.g.. lChabrie^l2003b . If 
the mass is decreased to 0.1 M© any derived velocity decreases 
by a factor of ~2 only and the characteristic distance increases 
by a factor of ~2. For a mass of 0.5 M©, the inferred char- 
acteristic infall distances are between 600 and 10^ AU, with 
Ser SMM4 and LI 527 representing the extreme cases. Three 
sources, IRAS15398, LI 157 and L1527, have such low infall 
velocities that the inferred infall radius is larger than the outer 
radius of the envelope as derived from the DUSTY modelling. 
Therefore, the entire envelope of these three sources is not in 
free fall, and the 10-K radius is adopted as the characteristic in- 
fall radius. Thus, the profiles trace large-scale infall ranging from 
several 100 AU to scales of > 1000 AU. 

Assuming spherical infall symmetry, the mass infall rates can 
be estimated through Minf = AnR^^^fimw f^inf where [i is the 



mean molecular weight (2.8 amu; iKauff'mann et"aDl2008h and 
^inf the density at the infall radius. For a free-falling envelope, 
this may be reformulated as 

M = 47T/dmiinM (IGM^^^f v^^^ 

= 1.3xlO-^M©yr-^f \( (1) 

\105cm-3;\o.5kms-V 

for a 0.5 M© enclosed mass (Menc ^ M^)- The density is taken 
from the relevant DUSTY model at the infall radius, Rm, and 
ranges from 7x10^ cm"^ (GSS30-IRS1) to 9x10^ cm"^ (Ser 
SMM4). Mass infall rates range from 3x10"^ M© yr'^ (GSS30- 
IRSl) to 2x10"^ M© yr"^ (LI 157). The very high mass infall rate 
inferred for LI 157 is a result of the very large infall radius (and 
therefore large surface area over which the infall takes place) 
and the very low velocity. As argued above, it is unlikely that 
IRAS15398, LI 157 and L1527 are in free fall, and the inferred 
properties are thus upper limits. We emphasize again, that the 
model used to derive these infall parameters is a toy model; to 
derive more accurate infall velocities and conditions it is neces- 
sary to run full radiative-transfer models of infalling envelopes 
which is beyond the scope of this paper. The values derived here 
serve primarily as order-of-magnitude estimates. 

The mass infall rate derived f or NGC1333-IRAS4A (9x 10"^ 
Mq yr"^) is very similar to what iDi Francesco et al. Idlooil) de- 
rived for the same source (1.1x10""^ M© yr"^). The infall radius 
is found to be nearly the same (2700 AU versus 3500 AU), and 
the density at that radius diff'ers by a factor of two (6.9x10^ cm"^ 
versus 3.0x10^ cm"^), thus accounting for the small diff'erence. 
For the case of NGC1333-IRAS4A, the interferometric observa- 
tions of H2CO and the H2O data presented here thus trace infall 
on the same physical scale, not limited by the optical depth of 
the H2O line. 

The location of the sources are marked in the correlation 
plots shown in Fig. [6l (with a square), but the sources do not 
stand out in terms of envelope or outflow parameters. Thus, the 
envelope parameters alone do not appear to determine whether 
a source shows an inverse P-Cygni profile or not, nor is it the 
particularly luminous or massive sources which show inverse P- 
Cygni profiles. Moreover, the geometry of the system does not 
seem to play a significant role in determining whether a source 
has an inverse P-Cygni profile; the LI 527 outflow is located 
close to the plane of the sky, that of NGC1333-IRAS4A is likely 
at an inclination of ~ 45° whereas the outflow from IRAS 15398 
is seen nearly pole-on. Further modelling of the sources showing 
inverse P-Cygni profiles is under way to properly constrain the 
infall parameters of the envelopes (Mottram et al. in prep.). 

4.1.5. Regular P-Cygni profiles 

Two sources have regular P-Cygni profiles superposed on the 
outflow component (IRAS12496 and L1551-IRS5; both Class 
I), with two more sources showing signs of P-Cygni profiles 
(L723 and HH46-IRS; the former being Class (0, the latter Class 
I; Fig. [TT]). A regular P-Cygni profile is a signature of expan- 
sion rather than infall. For the same reasons as above, the sig- 
natures are thought to be real and not the eff'ect of chopping 
into nearby emission or absorption or emission caused by a fore- 
ground cloud. 

A similar analysis as above can be performed to examine 
characteristic expansion velocities; the only change in the model 

^ Although L723 is formally classified as a Class object, it is also 
the Class object with the highest Tboi; therefore it may be close to 
transitioning to Class 1. 
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Table 5. Best-fit parameters of the infall modelling. 
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Notes. A free-fall collapse model does not describe these envelopes, and no values are derived for the characteristic infall radius, or mass infall 
rate. From Yildiz et al. (in prep.). Fixed parameter, see text. Mass infall rate derived through a sphere bounded by the 10-K radius. The 
infall rate is therefore a strict upper limit to the actual infall rate (see text). 

is that the infall velocity is negative. All profiles can be repro- 
duced by an expansion velocity of 0.5 km s"^ (Fig.[TT]), with the 
same uncertainty of ~ 0.2 km s"^ as found above for the infall 
velocities. A reversal of infall has been observed in the high- 
mass star-forming region SgrB2(M) close to the Galactic Centre, 
where it is proposed that the powerful outflows push out the sur- 
rounding e nvelope, while infa ll continues along certain preferred 
directions (lRolff'setal.ll2010l) . The high-mass star-forming re- 
gion W3-IRS5 shows the same phenomenon ('Chavarria et alj 
[20TQh and it is thus possible that this mechanism is at play on a 
smaller scale in these low-mass objects, assuming that the proto- 
stellar wind is strong enough to reverse infall in large parts of the 
envelope. In this respect, it is worth noting that the two sources 
showing the strongest P-Cygni profiles are both Class I sources 
where the envelope is less dense and easier to push aside. With a 
constant expansion velocity of 0.5 km s"^ the envelope will fully 
disperse on timescales of 10^ years, whi ch is not c onsistent with 
inferred Class I lifetimes of >10^ years (l Evans et al. 2009), indi- 
cating that the expansion is either a local phenomenon and does 
not involve the entire envelope, or that the phenomenon is tran- 
sient and turns on at a certain stage of the evolution when the 
outflow overcomes the infall motion of the envelope and then 
disperses the envelope in a short period. 

4.2. H2O versus CO 

One of the results from SWAS observations is that the H2O / 
CO abundance ratio increases with velocity, as shown for the Fig. 12. Ratio of H2O lio-loi / CO 3-2 line intensities as a func- 
case of Orion (Franklin et al. 2008). In our sample of 29 YSOs, tion of absolute velocity for all sources. All spectra have been 
the H2O lio-loi / CO 3-2 fine ratio increases with velocity, as centred at km s"^ The top panel shows the ratio for all line 
shown in Fig. [T2I from an initial value of -0.01 to -0.2. There is wings where both H2O and CO emission exceeds 4 cr. Red- 
no significant diff'erence between the line ratio measured in red- shifted line wings are shown in red and blue-shifted wings in 
shifted fine wings compared to those in blue-shifted wings. For blue. The average value is shown in black. The bottom panel 
L1448, where both H2O and CO emission are detected to i;>40 shows the average value with the standard deviation shown in 
km s"\ the line ratio remains constant at -0.3 for i;>20 km s"^ grey. The average values for the red- and blue-shifted line wings 
The n on-LTE statistical equil ibrium radiative transfer code ^re shown in red and blue, respectively. The average ratios are 
RADEX (Ivan der Tak e Tan i2007h is used to estimate the abun- t^^en m 1 km s bms. 
dance of H2O with respect to CO in the optically thin limit. The 
kinetic gas temperature of the H2O and CO emitting regions 

are assumed to be identical, and the derived abundance ratio is order-of-magnitude abundance ratios of H2O/CO. If the density 

shown for two diff'erent densities, ^h=10^ and 10^ cm"^ (Fig. is low, the observed range of line ratios corresponds to abun- 

\T3[ . Recent indep endent studies suggest that this assumption dance ratios of the order of 10"^ - 10"Mn the H2O and CO op- 

may not be valid (iBjerkeli et al.l l201l[ ISantangelo et al.l l2012l tically thin limit. For a standard abundance of CO of 10""^ with 

IVasta et aT]|2012h . but we use it here to illustrate the typical respect to H2, this translates to H2O abundances of 10"^-10"^ 
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Fig. 13. H2O / CO column density ratio as a function of kinetic 
gas temperature for three H2O optical depths, t=1, 10 and 25 
(CO is assumed to be optically thin in all cases). Contours are for 
intensity ratios of H2O lio-loi / CO 3-2. The top panel shows 
the ratio for a density of nu = 10^ cm"^ while the bottom panel 
is for a density of 10^ cm"^. 



relative to H2. For the high-density case, these numbers are a 
factor of 10 lower. 

The resulting H2O abundance is consistent with what is 
inferred elsewhere ( x(H20) ~ 10"^-1 0"^; Franklin et al. 2008^ 
iLefloch et al] l2010l iKristensen et al.l [2010b). Kri stensen et all 
(l201lh studied the L1448-MM outflow and analysed the H2O 
abundance in the line wings reaching the conclusion that the in- 
tegrated water emission is very optically thick (t>10) based on a 
tentative detection of broad H^^O emission. If the same exercise 
as above is performed but now for the optically thick case, higher 
H2O / CO abundance ratios are inferred for the high- velocity and 
high-density gas leading to H2O abundances of ~10"^ with re- 
spect to H2. Figure [13] shows the abundance ratio for three op- 
tical depths of the H2O line, r = 1, 10 and 25. The CO 3-2 
emission remains optically thin for all the considered physical 
conditions (tco^I)- Thus, to obtain even an approximate value 
for the H2O / CO abundance ratio, the optical depth of the H2O 
emission should always be taken into account. 

The weak H2O-CO intensity correlation (Fig.O is likely due 
to diff'erences in the density of the emitting material and the 
abundance ratio of H2O / CO. Because H2O is sub-thermally 
excited (the critical density for the lio-loi line is -10^ cm"^), 
the intensity scales as n(tl2)N(tl20) whereas the intensity of the 
thermally excited CO scales directly with N(CO) (see also Fig. 
[T3]) . Both of these proportionalities are only valid in the optically 
thin limit. The CO emissi on in the line w ings is usually close to 
optically thin (t<1; e.g. JCurtis et al.ll2Q10 . Yildiz et al. subm.), 
wher eas water is eff'ectively optically thin (Fig. [131 ISnell et aP 
I2OOO ). The H2O / CO abundance ratio is the same for all sources 
to within an order of magnitude (see above), and therefore the 
weak correlation is likely dominated by diff'erences in the den- 
sity. The density is estimated using the results from the DUSTY 
modelling (Appendix [C]). The correlation between H2O emis- 
sion and the density at 1000 AU is strong (r = 0.71; Fig. [6]) and 
does not degrade if the density is taken at 500 or 1500 AU. This 
further corroborates the result that variations in the H2O-CO cor- 
relation are dominated by density eff'ects. 

4.3. H2O and evolution: Class versus Class I 

A clear trend found in this study is the lower integrated 557 
GHz intensity of Class I sources compared to Class sources. 
Figure [14] compares the average Class spectrum with the av- 
erage Class I spectrum and there is a diff'erence of a factor of 8 
in integrated intensity between the two. For this comparison all 
spectra have been scaled to a common distance of 200 pc, and 
GSS30-IRS1, R CrA-IRS5 and HHIOO-IRS have not been in- 
cluded, since their excitation may not be representative of Class I 
sources (see Sect. 13. II) . Two evolutionary parameters have previ- 
ously been used in the litera ture: Tboi and L^-^^ / Menv (IChen et al.l 
ll995l[S^aceno et aPl 19961) . both parameters increasing in value 
with evolution. There is indeed an anti-correlation of H2O emis- 
sion with both evolutionary tracers (albeit weak, r = -0.60 and 
-0.61, respectively; Fig. [6]), thus showing that water emission 
decreases with evolutionary stage. 

Class I envelopes are generally characterised by lower 
masses and lower densities throughout, conclusions that are 
confirmed by the envelope modelling presented here. As the 
protostellar system evolves, the outflo w cavity opens, and th e 
outflow becomes less collimated (e.g.. lArce & Sargenil2006b . 
Furthermore, the material swept up by the outflow is less 
dense, because the envelope is less dense, and therefore the 
H2O intensity decreases. This has the efl'ect of decreasing the 
outflow force, as there is less material to interact with (e.g.. 
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Fig. 14. H2O spectra of Class and Class I sources averaged 
independently and compared. All spectra were shifted to km 
s"^ before averaging, rebinnedto 0.3 km s"^ channels and shifted 
to a common distance of 200 pc. GSS30-IRS1, R CrA-IRS5 and 
HHIOO-IRS are not included, see text. 

iBontemps et"aD I1996L iHogerheiide elaP Il997h . Therefore the 
decrease in H2O emission in Class I systems compared to Class 

systems may be directly linked to the weaker outflows, which 
in turn is linked to lower mass infall rates and envelope dispersal. 

It has previously been proposed that UV irradiation of the 
envelope may lead t o photodisso ciation of water, thereby lower- 
ing the abundance (Nisini et al. 2002). Als o, X-rays h ave been 
found to efl'ectively destroy water (Stau ber et al.l 1200 6). These 
mechanisms cannot be ruled out with the observations presented 
here, however no decrease in the H2O / CO ratio is observed in 
the Class I objects compared to the Class ones (see Fig. lA.ll in 
the Appendix). Because the bulk of H2O emission in the Class 

1 sources is associated with the medium component, it is only 
possible to analyse the H2O / CO ratio over a few km s"^ from 
the source velocity for most of the Class I outflows. To fully 
test this scenario, high S /N velocity-resolved observations of the 
photodissociation product OH are required, such as has recently 
been done for the high-mass object W3-IRS5 (iWampfler et al.l 
l20TTh . 

5. Conclusions 

Emission from the 0-H2O 1 lo-loi 557 GHz transition is detected 
in 28 low-mass Class and I YSOs from a sample of 29 sources 
within the WISH key programme on Herschel. Line profiles are 
complex and rich in diff'erent dynamical components. In partic- 
ular, the maximum width of the Class profiles extends beyond 
100 km s"^ Several dynamical components that have not been 
seen before in other molecular tracers show up in the H2O line 



profiles. These components include molecular "bullets" and reg- 
ular and inverse P-Cygni profiles, indicative of expansions and 
infall, respectively. In particular, the occurence of water bullets 
in more sources than was previously known indicates that this is 
a common phenomenon, and that water is one of the best tracers 
of these objects. 

Class and I sources appear to be quantitatively diff'er- 
ent in terms of water line profiles and their intensity. Class 
sources are associated with broad and medium components, both 
Gaussian in nature, whereas the Class I sources tend to be as- 
sociated more with the medium component. In sources where 
both a medium and broad component are detected, the medium 
component is spatially less extended than the broad component 
and located between us and the region causing the broad com- 
ponent. We speculate that the medium component is the ther- 
mal component of the 22 GHz H2O maser emission observed 
in these sources. If that is the case, the excitation conditions 
(r, n) are higher in the medium component than the broad. 
Infall is observed directly in six Class sources and one Class 
I, whereas expansion is observed primarily in Class I sources. 
Bullets appear exclusively in Class sources, but not in Class 
I sources. Absorption features are predominantly saturated in 
Class sources, but not in Class I sources, directly deomon- 
strating that the water abundance in the outer envelope is low, 
>10"^^-10"^. This confirms earlier ideas of the Class phase 
being dominated by prominent outflows, while at the same time 
being dominated by infall. By the time the YSO enters the Class 
I phase, outflow activity has decreased significantly, and infall 
has in some cases turned to large-scale expansion of the enve- 
lope. 

Modelling of the inverse P-Cygni profiles in Class objects 
leads to inferred infall velocities of ~ 0.1 - 1.0 km s"^ The pro- 
files are tracing large-scale infall (>1000 AU). Corresponding 
mass accretion rates are in the range of -10"^ - 10"^ M© yr"^ 
depending on source. The expanding Class I envelopes are mod- 
elled using the same prescription, leading to characteristic ex- 
pansion velocities of -0.5 km s"^ Expansion profiles have pre- 
viously been observed in high-mass star-forming regions, but 
this is the first time they are observed in several low-mass ob- 
jects. The dispersion timescale is short (-10"^ years) indicating 
that the expansion is either localized to small regions in the en- 
velope or that it is a phenomenon taking place at a certain stage 
of protostellar evolution when the outflow overcomes the infall 
motion of the envelope and then disperses the envelope in a short 
period. 

Water emission from the outflow is found to be controlled by 
the envelope density, which is used as a proxy for the density of 
the outflowing gas. Only small chemical difl'erences between the 
outflows are found, and any difl'erence is dominated by diff'er- 
ent excitation conditions. There is a tendency for more evolved 
Class I objects to appear weaker in H2O emission, an observa- 
tional fact that is linked to the lower density in the envelopes and 
the weaker outflow force. 

The water abundance with respect to CO is found to increase 
with velocity in all sources from 10"^ close to the source veloc- 
ity to 0.1 at velocities greater than 10 km s"^ from the source, 
assuming that the water and CO excitation conditions are simi- 
lar. The increase is uniform in both the red- and blue- shifted line 
wings. There is no systematic difl'erence between the H2O / CO 
abundance ratio in Class and I sources, although the Class I 
line wings are weaker and it is difficult to make this conclusion 
robust. 

The results presented here are one of the first steps to- 
wards an understanding of the role water plays in star formation. 
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Thanks to its high spatial resolution and sensitivity, HIFI is de- 
livering velocity-resolved water data for a much larger sample 
of lov^-mass protostars than has ever been done before on scales 
where the relevant physical processes occur. For the first time 
it is possible to study the evolution of water in star-forming re- 
gions, something that will continue to be done within the frame- 
work of the WISH programme. 
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10' 




Fig. A.l. Emission ratio of H2O lio-loi / CO 3-2 as a function 
of absolute velocity. All spectra have been centred at km s"^ 
The top panel shows the ratio for all line wings where both H2O 
and CO emission exceeds 4 cr. Line wings from Class sources 
are shown in red, wings from Class I sources in blue. The aver- 
age value is shown in black. The bottom panel shows the aver- 
age value with the standard deviation shown in grey. The aver- 
age values for the Class and Class I sources in red and blue, 
respectively. The average ratios are taken in 1 km s"^ bins. 

Appendix A: H2O supplementary material 

Table I A. II provides an overview of the specific observing mode 
used for each object, including (when applicable) the reference 
position used. Reference positions were primarily chosen from 
SWAS observations showing these regions to be emission free, 
and regions were in general checked against ground-based CO 
1-0 observations. 

Table IA.2I provides the decomposition of each profile into 
Gaussians. Only saturated absorption features are not listed here. 
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Table A.l. Reference positions chosen for each source and obs ids. 



Source 


Obs. mode" 


RA^ 


Ti 

Dec^ 


ARA^ 


ADec^ 


Obs id 


L1448-MM 


PSSW 






10' 


5' 


1342203202 


NGC1333-IRAS2A 


PSSW 


52? 3525 


31? 1665 






1342202067 


NGC1333-IRAS4A 


PSSW 


52?3525 


31? 1665 






1342202065 


NGC1333-IRAS4B 


PSSW 


52? 3525 


31? 1665 






1342202064 


L1527 


DBS 










1342192524 


CedllO-IRS4 


DBS 










1342201525 


BHR71 


PSSW 






3' 


0' 




IRAS 15398 


DBS 












L483 


PSSW 






y 


3' 


1342217691 


SerSMMl 


PSSW 


277? 4011 


1?3112 






1342208580 


Ser SMM4 


PSSW 


277? 4011 


1?3112 






1342208579 


Ser SMM3 


PSSW 


277? 4011 


1?3112 






1342208577 


L723 


DBS 










1342210077 


B335 


PSSW 


294? 2271 


7?4356 






1342196409 


L1157 


PSSW 






10' 


2' 




NGC1333-IRAS3 


PSSW 


52? 3525 


31? 1665 






1342202066 


L1489 


PSSW 


59?5100 


26? 6428 






1342203197 


L1551-IRS5 


PSSW 


68?0388 


17?9564 








TMRl 


DBS 










1342192525 


TMCIA 


DBS 










1342192527 


TMCl 


DBS 










1342192526 


HH46-IRS 


DBS 










1342196410 


IRAS 12496 


DBS 










1342201526 


GSS30-IRS1 


PSSW 






60' 


0' 


1342205302 


Ellas 29 


PSSW 






60' 


0' 


1342204011 


Oph-IRS63 


PSSW 






60' 


0' 


1342205304 


RN091 


PSSW 


249?5017 


-16?2703 






1342205297 


RCrA-IRS5A 


PSSW 


284?5913 


-36?5986 






1342215840 


HHIOO-IRS 


PSSW 


284?5913 


-36?5986 






1342215841 



Notes. Observations were done in position switch mode (PSSW) or dual-beam-switch mode (DBS). In the latter case, a throw of 3' was always 
used. Absolute coordinates (J2000) for the reference position. ^""^ Relative coordinates with respect to the source position. 
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Table A.2. Gaussian decomposition of the H2O line profiles. 



Source 


Component 




/ Tub dv 






'^source 






(K) 


(Kkms-^ 


(km s-^ 


(km s-i) 


(km s-^ 


L1448-MM 


Broad 


0.28 


15.68 


52.8 


14.5 


+4J 




EVH-R 


0.11 


2.17 


19.2 


61.9 






EHV-B 


0.07 


2.21 


30.8 


-47.0 






Narrow 


-0.19 


-0.99 


4.9 


4.8 




NGC1333-IRAS2A 


Broad 


0.12 


4.49 


36.4 


3.4 


+1.1 




Medium 


0.19 


1.61 


7.9 


9.6 




NGC1333-IRAS4A^ 


Broad 


0.43 


15.02 


33.0 


7.8 


+1.2 




Medium 


0.28 


3.00 


10.1 


0.0 




NGC1333-IRAS4B 


Broad 


0.41 


10.07 


23.6 


6.2 


+1.A 




Medium 


0.27 


1.43 


8.8 


5.0 




L1527 


Medium 


0.05 


0.99 


20.3 


5.6 


+5.9 


CedllO-IRS4 


Medium 


0.02 


0.35 


14.3 


4.9 


+3.5 


BHR71^ 


Broad 


0.15 


7.21 


44.6 


-2.7 


-4.4 




EVH-R 


0.05 


1.50 


29.0 


55.4 






EHV-B 


0.03 


0.82 


23.0 


-48.7 






Narrow 


0.04 


0.09 


1.8 


-12.3 




IRAS 15398^ 


Medium 


0.09 


1.78 


18.9 


0.6 


+5.1 


L483 


Medium 


0.09 


1.62 


15.3 


3.5 


+5.2 


Ser SMMl 


Broad 


0.16 


6.62 


39.0 


14.2 


+S.5 




Medium 


0.45 


6.46 


13.6 


7.4 




Ser SMM4^ 


Broad 


0.14 


7.72 


40.2 


0.7 


+S.0 


Ser SMM3 


Broad 


0.12 


5.92 


46.8 


3.2 


+1.6 




Medium 


-0.11 


-2.03 


17.3 


0.0 




L723'^ 


Broad 


0.04 


0.68 


17.4 


9.9 


+n.2 


B335 


Broad 


0.03 


1.15 


40.5 


6.7 


+^.4 




Narrow 


0.07 


0.08 


1.1 


9.3 




L1157^ 


Broad 


0.09 


2.64 


26.5 


4.2 


+2.6 




EVH-R 


0.01 


0.30 


20.5 


42.6 




NGC1333-IRAS3 


Broad 


0.19 


6.22 


30.1 


5.2 


+S.5 




Medium 


-0.14 


-0.90 


5.9 


4.9 




L1489 


Broad 


0.05 


0.87 


17.6 


4.1 


+1.2 




Narrow 


-0.06 


-0.07 


1.2 


5.8 




L1551-IRS5^ 


Broad 


0.02 


0.66 


28.1 


16.1 


+1.2 


TMRl 


Medium 


0.03 


0.36 


11.4 


4.4 


+6.3 




Narrow 


-0.04 


-0.06 


1.6 


5.7 




TMCIA 




<0.02 


<0.04 








TMCl 


Broad 


0.02 


0.49 


28.4 


3.4 


+5.2 


HH46-IRS^ 


Broad 


0.04 


1.23 


26.8 


9.8 


+5.2 


IRAS 12496^ 


Broad 


0.02 


0.58 


25.4 


1.8 


+2.3 


GSS30-IRS1 


Medium 


0.18 


3.32 


17.8 


1.9 


+2.S 


Elias29 


Medium 


0.13 


2.25 


16.7 


3.4 






Narrow 


-0.06 


-0.07 


1.1 


5.5 +5.0 




Oph-IRS63 


Medium 


0.01 


0.09 


6.6 


0.5 


+2.S 


RN091 


Narrow 


0.03 


0.50 


18.7 


6.8 


+5.0 


RCrA-IRS5A 


Medium 


0.12 


1.93 


15.2 


7.4 


+5.1 




Narrow 


0.36 


0.96 


2.5 


6.0 






Narrow 


-0.28 


-0.20 


0.7 


5.7 




HHIOO-IRS 


Medium 


0.05 


1.06 


19.1 


1.7 


+5.6 




Narrow 


0.20 


0.95 


4.5 


5.4 






Narrow 


0.35 


0.31 


0.8 


6.3 





Notes. All values are obtained from (multiple) Gaussian fits to each line profile. Non-saturated absorption components have negative values of 
and J Tmb dt;. Saturated absorption features are not included in this table. See Table [3] for an overview of what components are present 
in which source. Av corresponds to the FWHM of the component. Inverse P-Cygni source; component parameters not shown. Regular 
P-Cygni source; component parameters not shown. 
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Appendix B: CO supplementary material 



The CO 3-2 data are presented here. This includes the spec- 
tra overplotted on the H2O spectra (Figs. IB. II and IB. 21) . 
Magnifications to highlight the CO emission in the line wings 
are shown in Fig. IB. 31 where the spectra are shown on the same 
absolute scale as the H2O data in Figs. [T] and O The integrated 
intensities are provided in Table IB. ll 



Appendix C: DUSTY modelling 



All envelopes are modelled using the I D spherically symmet - 
ric dust radiative transfer code DUSTY (llvezic & Elitzur|[l997b . 
The procedure for m odelling the envelop es is t he same as out- 
lined previously in Schoier et~an (1200 2*) and 'J0rgense n et al.l 
(I2OO2I), this work being an extension and update of the latter. 
In this setup, the density profile is assumed to be a power-law 
with n oc r~P and the inner boundary, rin, of the envelope is set 
to be where the dust temperature is equal to 250 K. Three pa- 
rameters are varied: the slope of the density profile, p, the size 
of the envelope, Y = rout/^n, and the d ust optical depth at 100 
/ im, Ti oQ. The dust opacity is taken from lOssenkopf & HenningI 
(I1994I) . Table 1, column 5, the so-called 0H5, corresponding to 
dust grains with thin ice mantles. The radial extent and slope 
of emission as observed at 450 and 850 yum is modelled first, 
whereby p is determined. Secondly the spectral energy distribu- 
tion (SED) is fitted with the remaining free parameters, Y and 
Tioo. The envelope mass is directly related to tioo and is best 
constrained by the SCUBA 850 /im flux. 



sources with no available sub-mm continuum maps the emission 
was assumed to fall ofl' with a slope of 1.5, and the derived en- 
velope mass is therefore more uncertain due to poor sampling of 
the SED at long wavelengths (BHR71 and CedllO-IRS4, and 
HH46-IRS at 450 jim). Fixing the slope is most relevant for 
southern sources, for which little or no continuum data at 450 
and 850 yum exist. The SEDs were assembled from the litera- 
ture and consist of data points between 60 and 1300 jim. The 
long wavelength SED points (A > 200 jim) are obtained from 
SCUBA/LABOCA, Bolocam and SEST. The short-wavelength 
points now include the Spitzer-MIFS points not available previ- 
ously (lEvans et al and new Herschel-FACS data (Karska 
et al. in prep.. Green et al. in prep.). Model images were con- 
volved to the appropriate observing beam for comparison. To 
ensure that the mass of the envelo pe corresp onds to that mea- 
sured from the sub-mm fluxes re.g.. lShirlevet al. 2000), the 850 
and 450 jim fluxes were given high weight in the fit. 

The best-fit parameters along with the model results are 
listed in Table IC.ll An overview of the fitted data and results 
are shown in Figs. lC.lllC.3l S ources which have previously been 
analysed by J0rgensen etalj (^002) were re-analysed, and in 
general there is good agreement between the two sets of results. 
The uncertainty on p and tioo is of the order of 0.2-0.3, while it 
is somewhat higher for 7, of the order of -500. 

The envelope mass is measured either at the Tdust = 10 K 
radius or at the ^ = 10"^ cm"^ radius, depending on which is 
smaller. The envelope masses range from 0.04 M© (Elias 29) 
to ~16 Mq (Ser SMMl). These masses are in good agreement 
with masses derived from the sub-mm flux alone, e.g., through 
scaling of the 850 jim flux assuming a constant dust temperature 
( Shirley et al. 2000). The mass estimates presented here fall well 
within th e range of masses derived for these sources using other 



methods (Bontemps et al."l996', 'Hogerheiide et al 



IChen et a l. 2008, Enoch et al. 2009, J0rgensen etal 



199 
120091) . 



In this context it is important to note, that the modelling 
of the inner envelope is not realistic. Specifically, embedded 
disks are not taken into account and the model is expected to 
fail on scales smaller than a few 100 AU. In the Class I phase, 
disks can contr ibute significan tly to the sub-mm fluxes (e.g., 
|j0rgense n et al.ll2QQ7L ILommen e t al. 2008), and their contribu- 
tion has not been subtracted here. Therefore, the inferred enve- 
lope masses are upper limits, particularly for the Class I sources. 
The effect of an embedded disk is to add a point- source contri- 
bution to the flux at the center of the system and thus artificially 
increase the power-law density slope of the envelope. If the point 
source is subtracted, the /7-value would decrease as well as the 
overall envelope mass. Thus, the values of p and Menv presented 
here for Class I sources are overestimates. 



The emissio n profiles were ob tained from the SCUBA 
Legacy Archive (iDi Francesco et al.| [2008) or fro m LAB OCA 
(e.g., for the case of HH46-IRS: Ivan Kempen et aDl2009d) . For 
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Table B.l. Observed integrated CO 3-2 intensity in a W beam from the JCMT and APEX, and CO outflow force, Fco- 



Source 


rms"^ 




/ Tub dv 






Reference 




(mK) 


(K) 


(Kkms-^ 


(km s-i) 


(Mq yr"^ km s"^ 




L1448-MM 


0.07 


6.2 


120 


134 


9.6(-5) 


1,2,3 


NGC1333-IRAS2A 


0.11 


14.3 


106 


26 


2.4(-4) 


3 


NGC1333-IRAS4A 


0.06 


12.7 


69 


34 


2.3(-4) 


3 


NGC1333-IRAS4B 


0.04 


9.9 


46 


24 


1.2(-4) 


3 


L1527 


0.17 


12.4 


29 


7 


1.6(-4) 


1,4 


CedllO-IRS4 


0.77 


9.8 


50 


7 


9.6(-5) 


5 


BHR71 


0.05 


13.2 


85 


144 


2.6(-3) 


6 


IRAS 15398 


0.15 


8.8 


25 


10 


3.9(-5) 


5 


L483^ 


0.05 


7.7 


33 


17 


1.5(-5) 


1 


Ser SMMl 


0.10 


15.5 


133 


46 


l.l(-4) 


7 


Ser SMM4 


0.09 


25.8 


142 


43 


l.l(-4) 


7 


Ser SMM3 


0.09 


18.3 


191 


39 


3.8(-5) 


7 


L723 


0.07 


9.8 


64 


29 


6.0(-5) 


1 


B335^ 


0.38 


9.1 


44 


10 


2.4(-5) 


1 


L1157 


0.08 


8.0 


58 


43 


1.5(-4) 


8 


NGC1333-IRAS3 


0.09 


17.7 


116 


36 


6.2(-4) 


2,3 


L1489 


0.54 


5.7 


21 


5 


1.4(-5) 


4 


L1551-IRS5 


0.08 


8.2 


56 


26 


2.2(-3) 


2,4 


TMRl 


0.12 


6.2 


13 


7 


6.6(-6) 


1,4 


TMCIA 


0.13 


3.8 


35 


17 


8.7(-6) 


1,4 


TMCl 


0.14 


3.7 


19 


14 


2.6(-6) 


1,4 


HH46-IRS 


0.55 


19.0 


90 


17 


1.7(-4) 


5 


IRAS 12496 


0.42 


12.5 


85 


24 


e 




GSS30-IRS1 


0.10 


29.3 


120 


17 


e 




Elias29 


0.09 


19.8 


78 


17 


2.1(-6) 


1,9 


IRS63 


0.13 


11.5 


31 


12 


6.4(-7) 


1,9 


RN091 


0.09 


8.1 


23 


12 


4.2(-6) 


9 


RCrA-IRS5A 


0.35 


34.9 


174 


17 


e 




HHIOO-IRS 


0.39 


23.4 


120 


17 


4.4(-7) 


5 



Notes. Measured in 0.3 kms~^ bins. Line width measured at the 4-cr level. CO outflow force obtained from the literature. All values 
have been scaled to the distances used here. When multiple values were available, the average is used. Single-pointing data from the JCMT 
(beam-size is 15''). No values available in the literature. IRAS 12496 has an outflow seen nearly face-on, and so any value for Fco is uncertain 
at best; GSS30-IRS1 is located in a crowded region of Ophiuchus, and it is not possible to determine which outflow the source may be driving in 
single-dish CO 3-2 maps (van der Marel et al. in prep.). 

References. (1) Bontemps et al. (1996); (2) Cabrit & Bertout ( 1992); (3) Curtis et al. (201^; (4) lHogerheiide et all (119981) : (5) lvan Kempen et all 
(l2009h) : (6),BourkeetaL (.1997.) : m .Dionatos et al..(.201O) : f8) .Bachiller etal... (.2001) : (9) van der Marel et al. (in prep.). 
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Fig.B.l. Continuum- subtracted H2O lio-loi spectra of the observed Class sources as presented in Fig. [T]with CO 3-2 spectra 
scaled and overplotted in blue. The CO 3-2 scaling factor is written in blue for each spectrum. 
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Fig. B.2. Same as Fig.EUbut for Class I 
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Fig. B.3. CO 3-2 data for each source as obtained with the JCMT or APEX. The velocity and intensity scales are the same as in 
Figs.[I]and[2l 
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Table C.l. Best-fit DUSTY parameters and predictions. 



Fit 



Physical parameters 





n 
F 


Y 


'^lOO 


' in 


r(10 K) 




VI- 


f7(1000 ATI) 

1 lyLyjyjyj iiV^ j 


A/^(Ho) 














(M ) 




(rm~3) 
v^iii ) 


(rm~^) 
vCiii ; 


T 1 MM 


1 .J 


y\j\j 




90 7 
zu. / 




^ Q 
D.y 


l.^yy) 




7 7<'9^^ 




1 7 


son 




^S 




S 1 


'^.yyo ) 




J). / v-^-) ^ 


^^^-TR A<>14A 


1 R 
i .o 


1000 


7 7 






J.U 


^ UO) 
D. L\y ) 




1 0(94) 


IN VJ \^ 1 J J) J) -IJX/A.Ot-JJ 


1.4 


800 


'-T.D 


1 S 




^ 


9 0(Q) 


5 7(6) 


1 0(94) 


L1527 


0.9 


1200 


0.3 


5.4 




Q 


8.9(7) 


8 1(5) 


6 0(99) 


CedllO-IRS4 


1.4 


1400 


0.5 


4.1 


3 4(3) 


0.2 


8.5(8) 


3 9(5) 


1.2(23) 


BHR71 


1.7 


500 
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Fig. C.l. SEDs for each source (diamonds). The full line shows the best-fit SED. 
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Fig. C.2. 850 fim radial profile for each source (diamonds). The dashed line shows the beam profile, and the full line shows the 
best-fit radial profile. 




Fig. C.3. 450 fim radial profile for each source (diamonds). The dashed line shows the beam profile, and the full line shows the 
best-fit radial profile. 



